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ABSTRACT 

We present ~ 0"5 resolution near-infrared integral field spectroscopy of the Ha line emission of 14 z ~ 2 
UV-selected BM/BX galaxies, obtained with SINFONI at the ESO Very Large Telescope. The average Ha 
half-light radius is r-^n ~ 4^ kpc, and line emission is detected over > 20^,' kpc in several sources. In nine 
galaxies, we detect spatially-resolved velocity gradients, from 40 to 4 1 km s -1 over ~ 1 kpc. The observed 
kinematics of the larger systems are generally consistent with orbital motions. Four galaxies are well described 
by rotating disks with clumpy morphologies, and we extracted rotation curves out to radii > lOfc^ kpc. One 
or two galaxies exhibit signatures more consistent with mergers. Analyzing all 14 galaxies of the sample in 
the framework of rotating disks, we infer mean inclination- and beam-corrected maximum circular velocities 
of v c ~ 180 ±90 kms" 1 and dynamical masses from ~ 0.5 to 25 x 10 10 fc^o M Q within rw 2 . On average, the 
dynamical masses are consistent with estimates of photometric stellar masses assuming a Chabrier or Kroupa 
initial mass function (IMF) but too small for a 0.1 - 100 M Q Salpeter IMF. The specific angular momenta of 
our BM/BX galaxies are similar to those of local late-type galaxies. In the single disk framework, the specific 
angular momenta of the baryons are comparable to those of their dark matter halos. Extrapolating from the 
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average v c at IO/170 kpc, the virial mass of the typical halo of a galaxy in our sample is 10 
good agreement with previous estimates from the clustering and number density of the BM/BX population. 
Kinematic modeling of the three best cases implies a ratio of v c to local velocity dispersion v c /a ~ 2-4 
and accordingly a large geometric thickness. We argue that this suggests a mass accretion (alternatively, gas 
exhaustion) timescale of ~ 500 Myr. We also argue that if our BM/BX galaxies were initially gas rich, their 
clumpy disks will subsequently lose their angular momentum and form compact bulges on a timescale of 
~ 1 Gyr. In support of this scenario, the object with brightest A'-band magnitude and among the reddest in 
J — K colour in our sample (suggesting it is the most evolved) exhibits a maximum in [N II]/Ha ratio and a 
minimum in Ha equivalent width at the geometric/kinematic center, which we interpret as indicative of an 
inside-out metallicity and age gradient. 
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1. INTRODUCTION 

Over the last decade, deep surveys in the optical, infrared, 
and submillimeter wavebands have revealed a large diversity 
among high redshift galaxy populations. This reflects in part 
a diversity in properties such as mass, stellar populations, and 
gas and dust content. High redshift samples have become suf- 
ficiently complete to enable fairly robust estimates of the cos- 
mic evolution of the stellar mass and luminosity density, star 
formation history, and nuclear activity. About half the stellar 
mass density in the lo cal universe appears to have been as- 
semb l ed bv 7. ~ 1 (e.g.. | Dickinson et al.l2003tlRudnick et al.l 
2003; Fontan a et al. . 2003), after the major peak of star for- 
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mation and QSO activ ity at z ~ 2-2.5 (e.g. JFan et"aLll200U 
Cha pman etal. 120051) . 

It is, however, less clear exactly how galaxies were assem- 
bled. In particular, when and at what rate did galaxies of dif- 
ferent masses form? What is the connection between bulge 
and disk formation? Progress in solving these questions is 
hampered by our incomplete knowledge of various key phys- 
ical and observational aspects. For instance, we still have a 
poor understanding of the competing processes of cooling, 
angular momentum exchange and loss, and feedback from 
star formation and active galactic nuclei (AGN), which drive 
galaxy evolution. The links between the various high redshift 
populations, their origin, and their subsequent fates still have 
to be established. In the framework of semi-analytical models 
of structure formation, the description of the formation and 
evolution of galaxies remains uncertain because of the lack of 
observations constraining the parametrization of the complex 
physical processes involved. The most important ingredients 
include angular momentum, metallicity, stellar initial mass 
function (IMF), the spatial distribution of the gas and stars, 
feedback, and environmental effects such as merger frequency 
and harassment. A better knowledge of the spatially-resolved 
kinematics, morphologies, stellar populations, metallicities, 
and ionization state of the line-emitting gas within galaxies 
will greatly aid in addressing these issues. 
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In this context, we have begun a substantial program to 
investigate representative samples of galaxy populations at 
Z ~ 1 — 4 from integral field spectroscopy at near-infrared 
(NIR) wavelengths with SINFONI at the Very Large Tele- 
scope (VLT) of the European Southern Observatory (ESO): 
the high redshift galaxy Spectroscopic Imaging survey in the 
Near-infrared with SINFONI (SINS). For z = 1 -4, NIR ob- 
servations probe the rest-frame optical emission, which con- 
tains numerous diagnostic spectral features that have been 
well studied and calibrated from large samples at low redshift. 
SINFONI enables the full two-dimensional mapping of the 
spectral features, in seeing-limited or adaptive-optics mode. 
To obtain a panoramic view of the characteristic spatially- 
resolved properties, we target samples of galaxies selected 
by a variety of well-defined criteria. One of the class we 
have observed is the rest-frame UV-selected "BM/BX" galax- 
ies, which are efficiently selected in optical surveys from 
colour criteria analogous to those of the Lyman-break tech- 
nique but targeting the lower redshift rang e z ~ 1.5-2.5 
( Adelberge 'r et : al.ll2004t ISteidel et al. 1 120041) . In this paper, 
we present the results of a sample of 14 BM/BX objects, fo- 
cussing on the morphologies and kinematics traced by the Ha 
line emission. 

In § |3 we describe the observations and data reduction of 
our BM/BX galaxy sample. In § |5J we present the results 
of the Ha morphologies and kinematics for the entire sam- 
ple while in § |3 we present detailed case studies of three of 
the objects, which have the best quality two-dimensional in- 
formation. We discuss our results in §[5] and summarize the 
paper in §[6] Throughout, we assume a A-dominated cosmol- 
ogy with H = 70/! 70 kms" 1 Mpc" 1 , ^ m = 0.3, and A = 0.7, 
unless explicitely stated otherwise. For this cosmology, 1" 
corresponds to « 8.2 kpc at z = 2.2. Magnitudes are given in 
the Vega-based photometric system. 

2. OBSERVATIONS AND DATA REDUCTION 

We selected our BM/BX targets from the samples ob - 
served with lo ng- slit N IR spectroscopy of | Erb et al. I J 2003). 
Shanl evetal.l (E00l . ISteidel etal.l (120041) . and lErb et al. I 
(2006bJK Our criteria were a combination of target visibil- 
ity, night sky line avoidance for Ha, and Ha line flux, with 
emphasis on objects with spatially extended emission, spa- 
tially resolved velocity gradients, or high velocity dispersions 
based on the existing long-slit spectroscopy. To date, we have 
observed 14 objects: one BM source at z = 1.41 and 13 BX 
sources in the range z = 2.0 — 2.5, with a mean and median 
redshift of 2.24. Our sample includes one clearly identified 
pair, Q2346-BX404/405, with projected angular separation 
of 3. "6 (corresponding to 30 kpc at the redshift of the sys- 
tem). The optical and NIR phot ometric properties of our tar- 
gets are presented by Erb et al. (2006b) and are very similar 
to those of their - 100 BM/BX objects with long-slit NIR 
spectroscopy, spanning nearly the full range observed for the 
general BM/BX population in terms of 1Z and K s magnitudes 
as well as 1Z — K K and J — K s colours (see, e.g., Stei del et al.l 
2004; Reddv et al. 2005). For comparison with other surveys, 
three of the 1 1 targets of our sample with available /f s -band 
photometry have K s < 20 mag, and two of the eight with ad- 
ditional /-band photomet ry ha ve J-K s > 2.3 mag. We refer 
the reader to Er b et al.l ( 120 06b) for a more detailed discussion. 
Table[0lists the sample (along with the redshifts derived from 
Ha in §0 and summarizes the observations. 

We collected the data with SINFONI jE isenhaue ret al.l 
2003^ iBonnet etal.ll2004l Idillessen et al. 1120051) mounted at 



the Cassegrain focus of the VLT UT4 telescope during sev- 
eral observing campaigns between 2004 July and 2005 Oc- 
tober, as part of commissioning and guaranteed time observa - 
tions. SINFONI consists of SPIFFI ( Eise nhauer et al. I2003bl) . 
a NIR cryogenic integral field spectrometer equipped with a 
2K 2 HAWAII detector, an d a curvature-senso r adaptive optics 
(AO) module (MACAO; IBonnet et all 120031) . We used the 
H- or /f-band gratings to map the Ha line emission, depend- 
ing on the redshift of the galaxies. For Q2343-BX610, we 
complemented /f-band data with observations in the J and H 
bands to determine the abundance indicator "7^23" from the 
[O n] A 3727, [O III] AA4959,5007, and H/3 emission lines. 
We carried out most of the observations in seeing-limited 
mode with the 0"125 pixel" 1 scale, for which the total field of 
view (FOV) is 8" x 8". For two sources, the presence of a suf- 
ficiently bright nearby star enabled us to obtain AO-assisted 
observations. For Q1623-BX502 (guide star 8" away), we 
used the 0"05 pixel -1 scale with total FOV of 3. "2 x 3. "2. For 
Q1623-BX663, we employed the 0"125 pixel" 1 scale and 
locked on the guide star (13" away) for one-third of the to- 
tal integration time. In these configurations, SINFONI deliv- 
ers a nominal FWHM spectral resolution of R ~ 4500 in K 
(A = 1 .95 - 2.45 /xm), - 3000 in H (A = 1.45-1.85 /im), and 
~ 1900 in J (A = 1.10-1.40 ^m). 

We took the data in series of "AB" cycles, with typical 
nod throws of about half the SINFONI FOV so as to image 
the source in all frames, and jitter box widths of about one- 
tenth the FOV to minimize the number of redundant positions 
on the detector array. The individual exposure times varied 
from 300 s to 900 s depending mainly on the variability and 
intensity of the background and night sky line emission, in 
order to optimize the background subtraction and remain in 
the background-limited regime over the wavelength range of 
interest around the Ha line. The total on-source integration 
times ranged from 1 hr to 8 hr, driven by the surface bright- 
ness of the sources and our aim to map the kinematics out 
to large radii for those with the most extended Ha emission 
and the largest velocity gradients. We took exposures of stars 
(including the acquisition stars used for blind offsetting to the 
galaxies as well as late-B, A0V, and Gl V to G3V telluric stan- 
dards) interleaved with the target observations to monitor the 
seeing and the atmospheric transmission. 

The data reduction included the following steps. We per- 
formed wavelength calibration of the individual BM/BX ex- 
posures using arc lamp frames, which we verified and cor- 
rected, if necessary, using night sky lines present in the data. 
We subtracted the BM/BX frames pairwise (without spatial 
shifting) to remove the night sky line and background emis- 
sion. For some galaxies, we applied a more sophisticated 
sky subtraction procedure that involves scaling each transi- 
tion group of telluric OH lines separately in order to minimize 
residuals affecting the emission lines of interest of the source. 
We flat-fielded the BM/BX data with exposures of a halogen 
calibration lamp. We identified bad pixels based on dark and 
flat-field frames, and replaced them using interpolation. We 
then reconstructed (including a distortion correction) and reg- 
istered the three-dimensional data cubes, and averaged them 
together to produce the final reduced cube for each source. 
The data of the offset stars and telluric standards were re- 
duced in a similar way. Based on the reference stars data, 
the seeing was typically « 0"5 (FWHM) in the NIR through- 
out the observations. For Q1623-BX502, the use of AO at 
the 0."05 pixel" 1 scale improved the effective resolution to a 
FWHM of 0"3. From the night sky lines, we determined that 
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the effective spectral resolution of the reduced BM/BX data 
cubes corresponds to « 80 kms" 1 (FWHM). 

3. RESULTS 

We extracted maps of the velocity-integrated line flux, ve- 
locity field, and FWHM velocity width (hereafter Av) of Ha 
from the reduced data cubes of the BM/BX objects. We also 
extracted position-velocity (p-v) diagrams in a w 0."5-wide 
synthetic slit oriented along the kinematic major axis, identi- 
fied from the SINFONI data for each object. To increase the 
signal-to-noise (S/N) ratio, we smoothed some of the p-v dia- 
grams with a two-dimensional Gaussian of FWHM = 2 pixels, 
or 0"25 spatially and » 70 kms" 1 spectrally. For the ve- 
locity and velocity width maps, we smoothed the reduced 
cubes spatially with a Gaussian of FWHM = 3 pixels, result- 
ing in a typical effective spatial resolution of » 0"6 (0."34 for 
Q1623-BX502 observed in AO mode with the 0."05 pixel" 1 
scale). We determined the total line flux, the systemic red- 
shift, and the total Av from Ha by fitting a Gaussian profile to 
the spatially-integrated spectrum of each source. The veloc- 
ity fields and velocity FWHM maps were extracted from the 
smoothed data cubes by applying the same Gaussian fitting 
procedure to each spatial pixel as to the integrated spectrum 
9 . Except for Q2343-BX6 10 (which has a Zf-band magnitude 
of 19.2 mag), the continuum emission of the BM/BX objects 
in the sample is sufficiently weak so that no continuum sub- 
traction was necessary in making the linemaps. We measured 
v,, defined as half the maximum observed velocity gradient 
from Gaussian fits to spectra taken in apertures on opposite 
sides of each source along the kinematic major axis. We also 
estimated the intrinsic half-light radius r\ / 2 along the geomet- 
ric major axis from Gaussian fitting to the Ha linemaps. 

Table |3 lists the redshift, r\n, Av, and v r measurements. 
The half-light radii and FWHM line widths are intrinsic quan- 
tities, corrected for the effective spatial and spectral resolu- 
tion. Figure^shows the Ha linemaps and p-v diagrams along 
the kinematic major axis for all of the BM/BX galaxies in the 
sample. Figure[2]shows the two-dimensional spatial distribu- 
tions of the velocity and FWHM velocity width (here uncor- 
rected for the spectral resolution) of Ha for five of the galax- 
ies. Together with Q2343-BX610 (shown in Figure[3}, these 
are the "best cases" among the sample observed with SIN- 
FONI in terms of the combination of high S/N, large spatial 
extent and good spatial resolution, and large velocity gradient. 

3.1. Ha Spectra, Sizes, and Morphologies 

In most cases, the integrated spectra from the SINFONI 
data agre e well with the long-slit sp ectra published by 
lErbetal.l (120031 IShaplev et all (120041 and ISteidel etail 
J2004I) . These were obtained at lower spectral resolution 
(R ~ 1500) with NIRSPEC at the Keck II telescope except for 
one, which was obtained at comparable resolution (R ~ 3500) 
with ISAAC at VLT UT1. Differences (significant in a few 
cases only) can be attributed to difficulties in night sky line 
subtraction, S/N limitations for extended sources or broad 
lines, and slit orientation and positioning that missed part of 
the emission for the largest sources. In addition to Ha, we 
detected the [N II] A 6584 emission line in 12 of the 14 galax- 
ies. The integrated [N II] A6548/Ha line flux ratios (hereafter 
simply [N II]/Ha ratios) are listed in Table |2] (see also Fig- 
ureCJ. 

9 Throughout the paper, "velocity" and "velocity width" thus refer to the 
center and FWHM parameters of the Gaussians that best fit the observed line 
profiles. 



The Ha sizes and morphologies among the sample span a 
fairly wide range. The majority of the sources (1 1/14) are re- 
solved spatially. The intrinsic r^n along the major axis range 
from about 2 to 7 kpc, with a mean of 4.4 kpc (median of 
4.1 kpc). Three sources may be only marginally resolved, but 
the effective spatial PSF for their data is more uncertain. In 
the most extended sources, the Ha emission tends to be ir- 
regular and somewhat clumpy. In the six best cases (Figures|2] 
and[3}, the Ha emission is clearly detected over ~ 10-20 kpc. 
In all of them except Q1623 -BX528, the brightest Ha peak 
is located off the geometric center and the Ha morphology 
exhibits features reminiscent of a tadpole or chain-like shape. 

For two of the sources we observed, HS T broad-band op - 
tical imaging is available (see Figure 3 of lErb et al.ll2003l) . 
ACS and earlier WFPC2 data of SSA22a-MD41 reveal a 
very clumpy rest-frame UV morphology, with bright knots to 
the south-west and diffuse faint emission towards the north- 
east. The SINFONI linemap shows that the Ha emission fol- 
lows very well the overall distribution of the rest-frame UV 
light (accounting for smearing by the seeing in our ground- 
based data). Q1623-BX376a, in contrast, is one of the most 
compact in Ha among our sample. WFPC2 imaging indi- 
cates it is also compact in the rest-frame UV but it lies at the 
boundary between two detectors, making it difficult to assess 
its detailed spatial distribution. In these two cases for which 
we can make the comparison, the spatial distribution of the 
Ha and rest-frame UV emission is overall very similar. 

The large Ha sizes in several of the BM/BX sources ob- 
served with SINFONI may not be surprising in view of our 
selection criteria, which emphasized (but were not strictly 
restricted to) bright and/or extended targets as in dicated 
by previous long-slit spectroscopy (Krbet al.N2003L [2006b; 
Shaplev et al.1 120041 ISteidel et al.ll2004l) . Th e typical half- 
light radius of BX ga laxies is w 2 kpc (e.g. IShaplev et al. I 
2004; Erb et al. 1 120041) . It is perhaps more interesting that 
most of the extended sources in our sample do not appear 
regular o r centrally conce ntrated in Ha emission. From ACS 
imaging, lErb et al. lJ200l also found a prevalence of irregular 
and clumpy morphologies peaking off-center in the rest- frame 
UV for a sample of 13 BX galaxies, which were selected 
primarily to be elongated among their full spectroscopically- 
confirmed sample in the GOODS-North field. 

One may consid er these characteristics in ligh t of 
the results of lEl megreen. Elmepre en. & Hirstl ll2004al) and 
Elme green. Elm egreen. & Sheets ( 2004ljK These authors 
found from ACS imaging that the majority of faint blue 
galaxies at moderate to high redshift (z ~ 0.5 - 2) appear 
as "chains," "tadpoles," and "clum p-clusters" and exh ibit a 
lack o f central bulges (see also, e.g.. l Cowie. Hu. & Songaiial 
The analogy with the UV-selected BM/BX galaxies is 
instructive, although one must keep in mind the differences in 
selection and in redshift range in making this comparison. In 
addition, the rest-frame UV emission may not trace well the 
bulk of the stellar population and could miss a bulge compo- 
nent. In that respect, the recent study by Toft and coworkers 
may shed some new light (S. Toft et al. 2006, in preparation). 
These authors investigated the morphological variations be- 
tween the rest-frame UV and optical of z > 2 galaxies using 
deep HST ACS and NICMOS imaging at optical and NIR 
wavelengths. They found that the majority of their sources 
with blue rest-frame optical colours (as defined by observed 
NIR J-K < 2.3 colours) show little change in morphology, 
including those that are clumpy and irregular. In contrast, 
their red sources typically become more regular and often 
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bulge-dominated in the rest-frame optical. NIR /- a nd K-band 
photom etry exists for eight sources in our sample (E rb et all 
2006b) and all but two have J — K < 2.3. In view of this, 
and given that BM/ BX galaxies have in general J—K < 2.3 
colours (e.g. Reddv et al. 1 120051) . the clumpy and asymmet- 
ric morphologies in the SINFONI Ha linemaps and in HST 
optical imaging may reflect intrinsic features of the overall 
distribution of the stellar population with a lack of prominent 
central bulges. We return to this point in §@]and[5] 

3.2. Ha Kinematics 

From the integrated Ha line widths (corrected for instru- 
mental resolution), we derive a mean and rms velocity disper- 
sion of 130 ±50 kms" 1 for our BM/BX galaxies, with a large 
range from about 70 t o 240 kms" 1 . This is i n good agree- 
ment w ith the results of Er b et al.H2003U2004l)lShaplev et al.l 
(2004), and D. Erb et al. (2006, in preparation) from long-slit 
spectroscopy of various BM/BX samples. For nine galaxies, 
we find significant spatially-resolved velocity gradients (see 
Figure [Q. At first glance, this may not appear too surpris- 
ing. Indeed, four of the sources were preferentially chosen as 
having tilted Ha emission in long-slit spectra. In those cases, 
the new SINFONI data add two-dimensional information and 
higher spatial resolution. For the remaining five sources, how- 
ever, the existing long-slit data did not show any significant 
velocity gradient. 

For the six best spatially-resolved (since most extended) 
galaxies with two-dimensional kinematic information, the ve- 
locity fields appear smooth and ordered (see Figures|2]and|3}. 
Similar to local star-form ing galaxies (e.g., Elmouttie et al. 
H9981 Idordon et al.l 2000). a smooth kinematic structure — 
despite a clumpy and asymmetric Ha light distribution — 
probably is indicative of ordered orbital motions in a disk or a 
merger. With the exception of a few objects, our modest linear 
resolution (~ 4 kpc) compared to the source size does not al- 
low us to uniquely distinguish between orbital motions of two 
galaxies in a (widely separated) merger or rotation in a single 
disk. For the purpose of mass estimates using the virial the- 
orem discussed below, the difference between the merger and 
single disk cases, however, is unlikely to be very different. In 
analogy with the ultra-luminous infrared galaxies in the loca l 
universe (e.g., Downes & Solomon 1998; Genzel et al. 2001), 
the gas in an advanced merger probably lies mostly in a single 
plane given that it dissipates efficiently on order of a dynami- 
cal timescale and settles in a rotating disk- or ring-like struc- 
ture of moderately large scale at the center of the gravitational 
potential. For the same observed parameters and orientation 
of the orbital plane relative to the line of sight, a merger model 
yields a fact or < 2 times larger dynamic al mass t han a rotating 
disk model (Heisler, Tremaine, & Bahcall 1985). 

The fact that seven of the nine galaxies with spatially- 
resolved velocity structures exhibit the largest velocity gradi- 
ents approximately along the morphological major axis sup- 
ports a rotating disk i nterpr etation. This contrasts with the 
findings of Erb et al.l ll200 4) based on long-slit spectroscopy 
of mostly elongated BX sources in the GOODS-North field. 
The reason for this difference is unclear; it could possibly 
reflect in part a combination of selection biases for both 
samples and the comparatively lower velocity resolution (~ 
210 kms" 1 ) and angular resolution (mostly 0"7-0"9) for the 
lErb et al.l (120041) long-slit data. Stronger evidence for rotat- 
ing disks is that the velocity gradients along the major axis of 
SSA22a-MD41, Q2343-BX389, and Q2343-BX610 flat- 
ten at the largest radii probed by the SINFONI Ha linemaps, 



suggestive of the flat part of a disk rotation curve (see Fig- 
ure0J. These three sources also show a maximum in FWHM 
velocity width at their geometric and kinematic center (and 
not at the brightest Ha emission peak offset from the center; 
Figures|2]and[3}, which is a key property of rotating disks. Fi- 
nally, the velocity field of Q2343-BX610 appears to exhibit 
a two-dimensional "spider diagram" (Figure[5), a compelling 
signa ture of a classical rotating disk (Ivan der Kruit & Allenl 
1978). We discuss these three galaxies in more detail in §|4] 

In contrast, the reversal in velocities observed for Q1623- 
BX528 (Figure c ould be indicative o f a counter-rotating 
merger (e.g iMihos & Bothunl 119981 iTecza et al.l 120001 
iGenzelet al.l 1200 11) . The source Q1623-BX663 exhibits a 
smooth velocity field despite its complex surface brightness 
distribution, which could be either due to a relatively face- 
on disk or a merger. The largest FWHM velocity width (~ 
500 kms" 1 ) occurs off-center, towards the prominent north- 
eastern Ha peak. The latter characteristic is rather more con- 
sistent with the source being a merger remnant or disturbed 
spiral with two long tidal tails, similar t o the local analo g 
UGC 10214 (the "Tadpole Galaxy," e.g.. iTran etail 120031) . 
The larger velocity width towards the Ha peak, or "head," 
of Q1623-BX663 may then perhaps be attributed to an AGN 
located there, consistent with AGN spectral features present 
in its rest-frame UV spectrum (A. Shapley, private communi- 
cation). 

Keeping in mind that higher resolution observations will 
likely reveal a range of kinematic properties, which are per- 
haps very complex in any given source, the evidence for sev- 
eral of our best sources nevertheless encourages us to ana- 
lyze the entire sample in terms of a rotating disk toy model 
to derive first order estimates of circular velocities and dy- 
namical masses. We considered simple models consisting of 
azimuthally symmetric rotating disks, in which the input light 
distribution is assumed to trace the distribution of the mass 
surface density. Other input parameters include the total mass, 
the inclination angle i of the normal to the disk's plane with 
respect to the line of sight, and the vertical (z-axis) thickness 
of the disk from which the z component of the velocity disper- 
sion is computed under the assumption of hydrostatic equilib- 
rium in the limit of a thin and very extended disk. Further- 
more, the model allows for an additional one-dimensional, lo- 
cal isotropic velocity dispersion component <jq. The beam 
smearing due to seeing is taken into account by convolution 
of the inclined model with a two-dimensional Gaussian PSF 
of specified FWHM. 

We used these models to estimate the projected circular 
velocity at turnover, v c sin(/), from our measurements of v r 
(half the maximum observed velocity gradient) and Av (the 
FWHM velocity width) as follows. For rotating disk models 
that account for beam smearing and with a range of sizes and 
of local one-dimensional velocity dispersions appropriate for 
the BM/BX galaxies in our sample 10 , the ratio v c sin(z)/v r is 
~ 1.3. The Ha emission may preferentially trace the central 
parts of the galaxy inside the turnover radius, so that the ob- 
served velocity gradient could underestimate the intrinsic ve- 
locity difference at turnover. This "velocity gradient" method 
thus yields Vc dgrad sin(i) values that may represent lower limits 
to the actual v c sin(z). In the same framework and with similar 
ranges of parameters, v c sin(/)/Av is ~ 0.42. If there is any 
significant line broadening other than due to the projection 

10 For the purpose of this analysis, we assumed Gaussian light profiles but 
similar results are derived with exponential disks or ring-like distributions. 
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of the rotation and random motions related to the disk thick- 
ness (e.g., caused by outflows/winds or streaming motions), 
this "line width" method yields v" ldth sin(/) estimates that may 
constitute upper limits to the intrinsic v c sin(/). In agreement 
with these expectations, the velocity gradient method leads to 
values on average 0.8 times those obtained with the line width 
method. We find the best agreement for the six best resolved 
galaxies, for which the mean v ^ elgrad / v w idt h j s \ \ \y e adopted 
the average of the values from the two methods as final es- 
timates. For the 14 BM/BX galaxies studied here, v c sin(/) 
ranges from 60 to 250 tans" 1 . Assuming that the sources have 
random inclinations and correcting for 1 / (sin(0) = 1 .57 gives 
an average and rms (v c ) = 180 ±90 tans" 1 . Table|2]lists the 
dynamical masses enclosed within the Ha half-light radius, 
calculated as M& yn sin 2 (/) = [(v c sin(/)) 2 ?"i/2]/G, where G is 
the gravitational constant. For the three galaxies for which 
more detailed kinematic modeling is possible (§ |4}, the dy- 
namical masses estimated from the simple recipes above are 
close to (a factor of 0.82 on average) those derived from the 
modeling. 

3.3. Additional Considerations 

The new SINFONI data represent a significant step forward 
in our ability to analyze the kinematics of high redshift galax- 
ies because of the full two-dimensional spatial mapping at 
good angular and spectral resolution. A major advantage of 
integral field spectroscopy is that no a-priori assumption about 
the kinematic major axis needs to be made. Nevertheless, the 
interpretation of velocity profiles and maps obviously requires 
some caution for observations of distant systems where uncer- 
tain inclination and beam smearing due to seeing can play an 
important role (see, e.g., the discussion by E rb et al.l2 004). 

As stated above, we cannot rule out a merger inter- 
pretation for our BM/BX sources, except possibly for a 
few of them, mainly because of spatial resolution limita- 
tio ns. The simulations of in tegral field spectroscopy data 
bv lLaw. Steidel. & Erbl (12006) indicate that while it may still 
be difficult to distinguish between rotating disks and two- 
component mergers for faint and/or compact high redshift 
systems, such a distinction is however possible for extended 
sources if the data have sufficient angular and spectral res- 
olution, and sufficient S/N. This may be the case for the 
three best cases of our sample (discussed in § |4}, especially 
Q2343-BX610 for which the high quality data reveal key 
kinematic signatures favouring a rotating disk scenario. 

Other factors also need to be considered. For instance, 
outflows and winds may affect the Ha kinematics. In lo- 
cal starburst galaxies with "superwinds," the outflowing gas 
does not d ominate the emission line ki nematics along the ma- 
jor axis (|Lehnert& Heckman ll996al) . neither does it dom- 
inate the total high surfa ce brightness Ha line emission 
jLehnert & Heckmanll995lll996bS) . Even in loc al infrared lu- 
minou s starburst galaxies with large extinction, lArmus et al.l 
(1990) found that the extended emission line gas outside of 
a few kiloparsecs from the nucleus accounts for < 25% of 
the total Ha emission. For the luminous z = 2.56 dust-rich 
starburst galaxy SMMJ1401 1 +0252, our SINFONI data (see 
iTecza et al. 120041) indicate an outflow component but its con- 
tribution to the integrated Ha emission is < 10%. Differential 
extinction therefore does not seem in general to play a ma- 
jor role in the relative contribution of outflows compared to 
gas photoionized directly by the starburst. There is also lit- 
tle evidence for a significant amount of shock-heated gas in 



BM/BX gal axies (such a s high [N II]/Ha ratios, as discussed 
bv iLehn erj: <fe Heckman||1 996al see Table |2] and Er b et al.l 
120031 IShaolev et al. 120041 D. Erb et al. 2006, in preparation). 
It is thus unlikely that mechanically heated gas plays a signif- 
icant role in determining the kinematics along the major axis 
of BM/BX objects. 

4. CASE STUDIES: ROTATION CURVES OUT TO > 10 KPC AND 
METALLICITY AND STELLAR POPULATIONGRADIENTS 

For SS A22a- MD4 1 , Q2343 -BX389, and Q2343 - BX6 1 0, 
we carried out significantly deeper Ha observations than for 
the rest of the sample (integration times of 8, 5, and 6 hr, re- 
spectively). This has enabled us to study their kinematics on 
radial scales > 10 kpc and determine a detailed velocity and 
mass model out to the asymptotic part of the rotation curve. 
Figure|4]shows the velocity profile along the kinematic major 
axis derived from the SINFONI data of these three galaxies. 
For Q2343-BX610, the data further allowed us to investigate 
spatial variations of the metallicity and stellar population. 

We appl ied t he axisymmetric rotating disk models de- 
scribed in § !3.2l to fit the Ha kinematics of the galaxies. From 
the models, we then determined the following key properties: 
the inclination angle i, the total (r > 10 kpc) dynamical mass 
Mdyn and the dynamical mass within the half-light radius, the 
circular velocity at turnover radius v c and at half-light radius, 
the vertical scale height of the disk zu e , and the ratio v c /a 
at turnover radius (where a is the velocity dispersion in the 
z-direction). We explored a range of model parameters and 
fit the data by tuning the main parameters interactively. Our 
principal criterion for the goodness-of-fit was that the models 
reproduce the observed radial profile of the velocity along the 
major axis (Figure |4) as well as that of the FWHM velocity 
width. In addition, the models have to account for the major 
axis size (at FWHM or 20 percentile) and the minor to major 
axis aspect ratio. We otherwise gave little weight to the details 
of the spatial distribution of the Ha emission since it can be 
significantly affected, for instance, by extinction and by the 
actual distribution of H II region complexes within the galaxy. 
These, however, do not influence the overall kinematics traced 
by Ha. 

We explored exponential disk models as well as models 
where the light/mass distribution is in form of a ring, as is 
perhaps suggested by the observed Ha morphologies, which 
peak off-center. We found that both types of models can pro- 
vide good fits to the kinematic data. In terms of the surface 
brightness distributions, however, the exponential models are 
generally too centrally peaked and also have outer wings that 
are too large. Ring models do better in this respect but re- 
quire, in addition, a more extended component to reproduce 
the Ha emission at r > 1" that is observed in all three galax- 
ies. For the exponential disk models, a better fit is obtained 
if the disk has a central hole (of radius 0."1 -0."2), which is 
conceptually similar to a ring model. Exponential disks with 
scale lengths of 0."3 — 0."5, or rings with radii and thicknesses 
R w AR w 0"3-0."5 combined with a more extended com- 
ponent (e.g., a Gaussian of FWHM « 1"3), both match ade- 
quately the overall surface brightness distribution of the three 
galaxies. 

We note that of the 14 BM/BX galaxies we observed with 
SINFONI, only a very few (the most convincing case is 
Q2343-BX389) exhibit a classical "double-horn" integrated 
spectral profile, which is characteristic of motions in a thin 
azimuthally symmetric plane, as for rotating thin disks and 
some mergers. For such models to match the data satisfac- 
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torily, it is thus required that large random motions through- 
out the source wash out the double-horn pattern. In the case 
of rotating disks, another contributing factor may be a slowly 
rising rotation curve. The inevitable consequence of such ran- 
dom motions is that the planar structure must be geometrically 
thick. 

4.1. SSA22a-MD41 

SSA22a-MD41 is one of the largest sources among our 
BM/BX sample, with Ha emission detecte d ove r a projected 
diameter of 2f 9 (24 kpc). As discussed in § 13.11 the Ha light 
distribution follows well the rest-frame UV morphology and 
relative intensit ies seen i n deep H ST imaging. In the context 
of the lElmegre en et al.l J2004albl) analysis, SSA22a-MD41 
resembles a "tadpole" system with two bright knots dominat- 
ing one side of a clumpy disk or ring embedded in a more 
diffuse emission. There is no evidence for a central bulge. 

SSA22a-MD41 has the third largest projected velocity gra- 
dient in ou r sample, Vr ~ 14 kms" 1 , first detected in long- 
slit data bv lErb et al.l J2003I) . The observed maximum in ve- 
locity dispersion near the geometric and kinematic center of 
the source is w 110 kms" 1 . The two-dimensional velocity 
field is in excellent agreement with expectations for a rotating 
disk/ring model, with a smooth velocity gradient that is steep- 
est approximately along the morphological major axis (Fig- 
ure|2j. The most striking kinematic feature is the flattening of 
the velocity gradient at the largest radii where we still detect 
Ha emission, r rj l."5 (« 12 kpc) from the geometric center 
(Figure^. This feature is consistent with the asymptotically 
flat part of a classical rotation curve. 

The Ha kinematics and light distribution are very well fit by 
the rotating disk models described above. The total dynamical 
mass inferred (within the radius of 12 kpc probed by our data) 
is M dyn = (6.2 ± 1.5) x 10 10 M Q , where the largest source of 
uncertainty comes from the choice of light distribution (expo- 
nential disk with central hole or ring) and from the uncertainty 
in inclination angle (i = 71° ±5°). As would be expected 
from their larger spatial wings, exponential disk models give 
masses near the upper end of this range, while ring models 
give masses near the lower end. About 60% of this total mass 
is enclosed within the Ha half-light radius of rj/ 2 ~ 6.0 kpc. 
The intrinsic (i.e., model) rotation curve rises rapidly to a 
maximum of rs 165 kms" 1 , near rm, with a slow fall-off at 
larger radii. In the framework of our simple constant (expo- 
nential) z scale height plus constant additional isotropic ve- 
locity dispersion models, the best fits require a z scale height 
°f Zi/e ~ 0."17 and a ratio of peak rotational velocity to total 
z-velocity dispersion of v c /a = 2.2 ±0.5. Clearly the disk is 
geometrically thick and highly turbulent. 

4.2. Q2343-BX389 

Q2343 -BX389 has the largest projected velocity gradient 
among our sample, v r rs 235 kms" 1 , and Ha is detected over 
a projected diameter of 2f7 (22 kpc). The observed peak in 
velocity dispersion near the geometric and kinematic center 
is very large, rj 180 kms" 1 . The Ha m orphology is remi- 
niscent of "chain" galaxies (Elme green et al.ll2004albl) . with 
three main knots forming a linear structure of half-light radius 
r\n ~ 7 kpc. Both the two-dimensional velocity field and the 
radial velocity profile are as expected for a nearly edge-on ro- 
tating disk (Figures |2]and|4j. The velocity gradient is steepest 
along the morphological major axis and flattens at the largest 
radii probed by the data, r w 1 "2 (w 10 kpc). 



The kinematic features and the overall light distribution 
of Ha are well reproduced by disk models with an inclina- 
tion angle of i = 80° ±5°. The total dynamical mass derived 
from the models is M dyn = (1.95 ±0.5) x 10 11 M Q , of which 
about 70% is enclosed within r\h = 7.0 kpc. Again, uncer- 
tainties are dominated by those of the light distribution and 
inclination but in this case, the difference between exponen- 
tial disk and ring models is less than for SSA22a- MD41. 
Q2343-BX389 has the largest derived dynamical mass of our 
sample. The intrinsic rotation curve reaches a maximum of 
ri 300 kms" 1 , near the half-light radius, and decreases slowly 
at larger radii. Good fits also require the combination of a 
substantial z scale height (zi/ e ~ Of 14) and corresponding z- 
velocity dispersion, as well as an additional isotropic velocity 
dispersion of 85 kms" 1 . This implies v c /a = 3.3 ±0.7. 

4.3. Q2343-BX610 

Q2343-BX610 has the second largest projected size and 
velocity gradient among our sample, with Ha detected over 
2f4 (20 kpc) and v r « 165 kms" 1 . The half-light radius is 
r\j2 ~5.4 kpc. The observed peak in velocity dispersion 
near the geometric and kinematic center is Ri 160 kms" 1 . 
Q2343-BX610 is remarkably similar to SSA22a-MD41 in 
its morphological ("tadpole" shape) and kinematic features. 
Its observed rotation curve also appears to flatten at radii > 1" 
(> 8 kpc) from the geometric center (Figure^. This galaxy 
constitutes our best data set in terms of S/N and kinematic 
detail. This enabled us to map the two-dimensional distribu- 
tion of velocity and FWHM velocity width with good signifi- 
cance over a large fraction of the galaxy, encompassing about 
2 square arcsec (or 16 spatial resolution elements). 

The two-dimensional v and Av distributions are in excellent 
agreement with expectations for a moderately inclined rotat- 
ing disk. The velocity map shows the telltale "spider diagram" 
pattern (Figure[5} and the velocity dispersion clearly peaks at 
the geometric center of the source (FigureEJ, both basic char- 
acteristics of rotating disks (van der Kruit & Allen 1978). For 
this galaxy, we also detected the continuum in the /T-band, 
tracing the rest-frame optical continuum emission from stars 
that likely dominate the mass in the central regions. The spa- 
tial coincidence of the continuum emission peak and of the 
kinematic center from the velocity map lends further support 
for the rotating disk scenario. 

Rotating disk models (either exponential disks with a cen- 
tral hole or rings) provide excellent fits to the data of Q2343 - 
BX610, including in this case the two-dimensional velocity 
field of Figure |3 The axis ratio of the emission indicates an 
inclination angle of i = 63° ±5°. The total dynamical mass 
of this disk is Md yri = (9.2 ± 1.4) x 10 10 M , of which 75% 
is enclosed within the half-light radius. The difference be- 
tween masses from exponential disks and from ring models 
is small. The peak rotation velocity is v c = 225 kms" 1 , near 
the half-light radius. As for the other two sources, our sim- 
ple models require a substantial z scale height (z\i e ~ 0."16) 
and associated z velocity dispersion, and an additional large 
isotropic velocity dispersion (77 kms" 1 ). The inferred v c /a 
ratio is 2.4 ±0.6. 

4.4. The Dynamical Center, K-band Surface Brightness, and 
Metallicity Gradient in Q2343-BX610 

Q2343 -BX6 1 has the brightest K-hanA magnitude and the 
highest [N II]/Ha ratio of the BM/BX sample studied here, 
with total K = 19.2 mag (D. Erb et al. 2006, in preparation) 
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and [N II]/Ha = 0.40 (Table|2j. Its brightness has allowed us 
to investigate the relative spatial distribution of the line and 
continuum emission, and their relation to the emission line 
kinematics. 

Figure [5] shows the detailed properties of Q2343-BX610, 
including maps of the velocity-integrated Ha and [N II] line 
flux, the [N II]/Ha ratio, and the line-free Zf -band continuum 
tracing the rest-frame optical continuum. The [N II]/Ha ratio 
clearly varies spatially, reaching a maximum near the geomet- 
ric center of the continuum emission. The peaks of [N II]/Ha 
ratio and ZT-band continuum emission are within ~ 0"2 of 
each other. While there is Ha emission detected at this po- 
sition, there clearly is no local maximum. The [N II]/Ha 
ratio measured in a 0."5-diameter circular aperture centered 
where the ratio peaks is 0.55 ±0.12. Within an annulus of 
inner diameter 0."5 and outer diameter l."5, we estimate an 
average ratio of 0.38 ± 0.04 for the rest of the galaxy. The 
maximum value of [N II]/Ha is rather high but similar ratios 
have bee n observed in H II regions of nearby spiral galax- 
ies (e.g.. lRvderf ll995). Furthermore, there is no evidence 
for the presence of an AGN in a deep optical spectrum ob- 
tained with LRIS at the Keck I telescope (D. Erb et al. 2006, in 
preparation). In support of this result, our SINFONI //-band 
data give an integrated ratio of [O III] A5007/H/3 = 1 .5 ± 0.4, 
which, combined with the integrated [N II]/Ha ratio, places 
Q2343-BX610 in a region still populated by non-AGN 
galaxies and consistent with normal star-forming galaxies in 
the [O III] A 5007/H/3 versus \N II] A6584/Ha diagnostic di- 
agram (e.g.. iLehnert & Heckmanlll994l iKewlev et al. 1120011 
Kauffman n et al.ll2003l) . We conclude that the spatial varia- 
tions of the [N II]/Ha ratio and continuum emission are most 
likely due to gradients in metallicity and stellar populations. 

The relationship bet ween oxygen abundan ce and [N II]/Ha 
flux ratio proposed by Pettini & Pagel (2004, "N2" index), 

12 + log(0/H) = 8.90 + 0.57 x log([NII]/Ha), (1) 

implies 12 + log(0/H) = 8.75 ±0.04 at the location of the 
[N II]/Ha peak and 8.66 ±0.03 for the rest of the galaxy. 
For an assumed sol ar value of 12 + log(0/H) = 8.65 (e.g. 
Asplu ncfet al.l l2004). we infer an oxygen abundance about 
1.25 times solar towards the center of Q2343-BX610 and 
approximately solar for the rest of the galaxy. Since the [N II] 
emission, and thus the [N II]/ Hq ratio, tends to satur ate at so- 
lar metallicity and above (see Pettini & Pagel 2004, and ref- 
erences therein), the abundance gradient in Q2343-BX610 
could be even higher than implied by the application of equa- 
tion 0. We also used our additional J- and //-band data 
for this source to estima te the globa l abundance by apply- 
ing the R 2 3 method (Page l et al. 119791) as we ll as the "03N2" 
index calibration of Pettini & Pagel (2004), with appropri- 
ate corrections for e xtinction and ionization effects (as in, 
e.g., iTecza et al"ll2004l) . Integrated over the source, the R23 
method implies 12 + log(0/H) = 8.54 and the 03N2 index 
yields 12+log(0/H) = 8.55, both of which are consistent with 
the estimate based on the N2 index 1 1 . The observed velocity 
of the higher metallicity region indicates that it is located near 
the kinematic center of the galaxy. In addition, the relative 
spatial distributions of the ZT-band and Ha emission imply 
a significant decrease in the Ha equivalent width (EW) to- 
wards the center. The rest-frame EWs within a 0"5-diameter 
aperture centered on the ZT-band peak and integrated over the 

1 1 The S/N ratio of the J- and //-band data prevents us from measuring 
reliably spatial gradients in line ratios and metallicity indicators. 



galaxy are = 37 ± 3 A and 53 ± 2 A, respectively, sugges- 
tive of an increase in the age of the stellar population towards 
the ZT-band continuum peak (see, e.g., the m odel predictions 
from STARBURST99; Leither eTetaLlT 999). 

The SINFONI data of Q2343-BX610 thus indicate that 
the Zf-band continuum peaks very close to the kinematic cen- 
ter of the galaxy, while the Ha emission, tracing intense 
star-forming sites, peaks off-center. The observed gradi- 
ents in [N II]/Ha ratio and Ha EW suggest significant spa- 
tial variations in metallicity and stellar ages, with higher 
metallicities and older ages towards the dynamical center. 
Overall, Q2343-BX610 appears to be one of the most 
evolved galaxies among our sample. In addition to its high 
[N II]/Ha and bright ZT-band flux, it has a red J-K = 2.24 
colour, very close to the criterion for selecting z > 2 "Dis- 
tant Red Galaxies" (DRGs, w ith J-K > 2.3; Fra nx et all 
2003; van Dokkum et al. 2003), which are found to be gen- 
erally more evolved than UV-selected galaxies at similar red- 
shifts (iForster Schreiber et al. ll2004ULabbe et al.ll2005l) . The 
overlap between BM/BX galaxies and DRGs is small at K < 
21 mag (Reddv et al. 2005), but the metallicity estimates for 
Q2343-BX61 are consistent w ith the properties both of 
DRGs (Ivan Dokkum et al~l 120041) as well as of the brighter 
part o f the BM/BX population {K < 20.5; IShaplev et all 
2004). 

4.5. Mass-Metallicity Relation 

IShaplev et al. ] d200l discussed the first evidence for 
(super-)solar metallicities and the possible existence of a 
luminosity(mass)-metallicity relationship for BM/BX galax- 
ies. The measurements of the [ N II]/ Ha ratio and the dynami- 
cal mass estimates inferred in § !3.2l for our 14 BM/BX galax- 
ies (see Table |3 show considerable scatter, although the two 
objects with the largest dynamical masses have higher than 
average [N II]/Ha ratio. Converting the [N Ill /Ha ratios to 
O/H abundances with the Pettini & Pagel (2004) relation (see 
e q- IDX we can compare our O/H and Md yn estimates with 
the local relationship bet ween m etallicity and stellar mass de- 
rived by Tremonti et al. (2004) using the Sloan Digital Sky 
Survey. Here, it must be borne in mind that a scaling of our 
dynamical masses within the half-light radius r\ n to total dy- 
namical masses will lead to an upward correction, and from 
dynamical to stellar masses will lead to an (unkown) down- 
ward correction in the mass coordinate. Nevertheless, the data 
suggest that, overall, our BM/BX galaxies lie at lower metal- 
licity (by roughly 0.4 dex) for a given mass with respect to the 
local re lationship. This is consistent with the earlie r results of 
Shanl ev et al. U20041) as well as the recent work bv lErb et al. I 
( 2006a), who discuss in detail the stellar mass-metallicity re- 
lation at z ~ 2 based on a large BM/BX sample (including 
the majority of our SINFONI targets). The scatter and un- 
certainties of our [N II]/Ha ratios, the small mass range cov- 
ered by our sample, and possible systematic differences be- 
tween abundances derived from different methods (see, e.g., 
iTremonti et alll2004l IShaplev et al1l200H lErb et al.l l2006a. 
and references therein) make it difficult to assess the presence 
of a mass-metallicity trend from our data. 

5. DISCUSSION 

While the spatial resolution of our ground-based seeing- 
limited SINFONI data is still somewhat coarse, they provide 
us with powerful new tools to investigate the properties and 
possible evolution of high-redshift galaxies based on their 
kinematics. In this section, we explore some implications of 
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our results. For simplicity, we exclude the z ~ 1 .41 BM object 
of our sample and focus on the 13 z = 2-2.5 BX objects in 
the following analysis. 

5.1. Angular Momentum 

Within the context of the ACDM model, the large veloc- 
ity gradients within ~ 10 kpc we observe for the BX galax- 
ies of the SINFONI sample are derived from the tidal torques 
gener ated by merging dark matter halos ( Peebles 1969i lWhitel 
1984). If this hypothesis is correct, then these gradients are 
directly linked to the spin and specific angular momentum of 
the dark matter halo. Moreover, the amplitude of the gradients 
over large radii imply that the galaxies in our sample have sig- 
nificant amounts of specific angular momentum. Combining 
the inferred sizes and rotation velocities, we can make sim- 
ple estimates of the specific angular momenta of the baryons. 
An estimate of the angular momentum depends on the mass- 
weighted velocity distribution of the galaxies for which we 
have only limited constraints. Assuming a simple model, 
the angular momentum of the baryons can be estimated as 
j = (3 r l / 2 v c , where (3 is a dimensionless number that depends 
on the geometry and spatial distribution of the mass. For a 
disk-like geometry, values of (3 range between 0.7 for a con- 
stant density disk with outer radius r\n and 2 for an expo- 
nential disk of scale length r\n. Adopting (3=1, the average 
specific angular momentum for the BX sources of our sam- 
ple is then j^ x ~ 900 kms" 1 kpc. The larger rotation veloci- 
ties and sizes inferred for SSA22a-MD41, Q2343-BX610, 
and Q2343-BX389 imply higher specific angular momenta 
of j% x ~ 1000-2000 kms" 1 kpc. Interestingly, these val- 
ues are com parable to those of local late-type spirals (e.g., 
Aba di et al. 1200 3. and references therein). 

The theoretical explanation for local disk galaxies having 
significant angular momenta compared to that calculated for 
their dark matter halos is that there has been little dissipa- 
tion and net transfer of angular momentu m between the disk' s 
baryons and the dark matter halo (e.g., lAbadi et al.ll2003h . 
How does the specific angular momentum that we estimate 
for our sample of BX galaxies compare with that expected for 
dark matter halos at z ~ 2? In the c ontext of the discu ssion by 
IMo. Mao. & White! (11998! see also lWhite & FrenUl991l) . the 
specific angular momentum of the halos for our BX sample 
(with (z) = 2.2) is roughly given by: 

; h BX « 10 2 - 8 A005 f A* ( —) 1 kms" 1 kpc, 

Jh 0.05 V 180 kms" 1 / V 3-2 7 07 P 

(2) 

where A0.05 = A/0.05 is the normalized spin parameter of the 
dark matter halo and where we have assumed A = 0.05 since 
this is approximately the most freq uent value in th e spin pa- 
rameter probability distribution (see lMo et al. 1 19981) . and v c is 
the halo circular velocity at virializatio n rh. v i r in the classical 
spher ical collapse model (IWhite & Frenkl l 991tlMo & White! 
2002). Due to the variations in the spin parameter and to pos- 
sible differences in the mass distribution of individual halos, 
at constant halo circular velocity the angular mome ntum of a 
dark m atter halo can vary by a factor of a few (see Mo et al.l 
1998, and references therein). This large range and the uncer- 
tainty of our estimates of the specific angular momenta means 
that we can only make an approximate comparison with the- 
ory. If we hypothesize that the circular velocities at turnover 
of our BX galaxies are repr esent ative of the halo virial circular 
velocities (as we argue in § I5.2> . then these velocities and the 
range of estimated j® x values imply j® x / j^ x ~ 1. Estimating 



this ratio to be of order of unity implies that it is indeed likely 
that the baryons at small radii did acquire their specific angu- 
lar momentum during the collapse of their parent dark matter 
halos. Given the uncertainties, and more speculatively, this 
would suggest that the baryons do not transfer their angular 
momentum to the dark matter upon collapsing from the virial 
radius to the disk scale length of a few kiloparsecs. 

5.2. Constraints on the Properties of Dark Matter Halos 

Given that the specific angular momentum of the dark mat- 
ter halo appears to be comparable to that of the baryonic disk 
for the BX galaxies in the sample studied here, it is, in prin- 
ciple, possible to constrain the properties of their dark matter 
halos using the observed kinematics. 

In their analysis of the spatial distribution of UV-selected 
galaxies at 1 .4 < z < 3.5 in 21 disjoin t fields with total area 
of 0.81 deg 2 , Adelberge retal.l (l2b05a) deduced a correlation 
length for BX objects with (z) = 2.2 of 4.2 ± 0.5 hr 1 Mpc. 
Through comparison with ACDM num erical simulations of 
large-scale structure, Adelbe rger et al.l then inferred that the 
estimated correlation length corres ponds to halo masses of 
M? x ~ 10 11 8 -lO 122 M . Similarlv. lAdelberger et al.ld2004l 
2005a) found that the observed comoving number density 
of BX objects is consistent with the estimated clustering 
strength, meaning there are sufficient halos to host galaxies 
of these halo masses. 

If the circular velocity determined on < 12 kpc scales from 
our BX data is a reasonable estimate of the circular velocity 
of the halo the galaxy resides in, we can directly estimate the 
halo mass using the basic relationship between halo mass Mj, 
and circular velocity at virializatio n v c .(rh v \r) in the classical 
spher ical collapse model JWhite & Frenk 199ll lMo & White! 
2002). For our BX sample, this yields an average andrms 

For SSA22a-MD41, Q2343-BX610, and Q2343-BX389, 
the average Mj, is about a factor of two higher than for the 
whole sample. These Mh estimates based on the kinematics 
are remarkabl y close to those obtained by Adelberge r et al.l 
(2004, 2005a) from the clustering properties and comoving 
number density of the BX population. This supports our hy- 
pothesis that the circular velocities we derived for our BX 
sample may indeed be simila r to those at the halo virial ra- 
dius. T he analysis of lWhit e & Frenkl ill 9911) and lMo & White! 
(2002) also implies a virial radius for the halos of our BX 
galaxies of 

For SSA22a-MD41, Q2343-BX610, and Q2343-BX389, 
the virial radii are about 30% larger, or ~ 100 kpc. For these 
three cases, for which our data show a flattening of the rota- 
tion curve at r max ~ 10 kpc, we thus i nfer r max / n, V \ T y 0. 1 . 
At a radius of 0. 1 rh, v h, simulations (e.g. lHavashi et al. 12004) 
show that v c (r)/v c (rh,vir) ~ 0.9, supporting our contention that 
the circular velocities are a reasonable representation of the 
halo virial velocity for our BX objects. Since these simula- 
tions do not include the condensation of baryons at the center 
of these halos — which obviously happened for the galaxies 
studied here — the values of v c (r) inferred from these models 
constitute lower limits. Overall, the simple analysis above im- 
plies that the circular velocities we derived for our BX sample 



Integral field spectroscopy of z ~ 2 galaxies 



9 



are very likely good estimates of the true circular velocity of 
their halos. 

We note that in a second paper on the spatial clustering fo- 
cussing on z ~ 2 B X galaxies in four fields with NIR imaging, 
Adelberge r et al.l J20 05b) found a substantially higher corre- 
lation length for those brighter than K = 20.5 mag compared 
to the overall population (by a factor of about 2.4). However, 
as the authors point out, much of the signal leading to this high 
correlation length comes from a single field, whi ch may bias 
the av erage result (see also the discussion by Steid elet al.l 
120051 who found a significant redshift ov erdensity for BX 
source s in another of the four fields used by Adelberge r et al.l 
2005B. 

The estimates of the halo masses and galaxy dynam- 
ical masses of our BX sample imply that only a frac- 
tion of the baryons are associated with the galaxies (see 
lAdelberger et al.l2005al for a related discussion). By combin- 
ing our derived dark matter halo masses with the dynamical 
masses, we can estimate the amount of the possible total bary- 
onic matter in place in our BX galaxies. For SSA22a-MD41, 
Q2343-BX610, and Q2343-BX389, our dynamical model- 
ing (§14} leads toM dyn (r max ) w (0.6-2) x 10 11 M©. The halo 
masses derived above assuming virial equilibrium are about a 
factor of 10 higher, i.e., the ratio M dyn /Mh w 0.1. For com- 
parison, the currently best value for the fracti on of bary onic to 
total mass is S7b / Sl m ~ 0.17 (e.g. JSpergel et al.120 03). Taken 
literally, our results would suggest that even if the dark mat- 
ter content is negligible inside a radius of r max , only about 
half or less of the baryons associated with the dark matter 
halo are actually within r max . For the Milky Way, the frac- 
tion of total baryonic to d ark matter mass is 0.08 ±0.01 (e.g., 
ICardone & Serenoll2005l) . The similarity between the Milky 
Way and the three BX galaxies considered here suggest that 
at epochs as early as z ~ 2, at least some of the BX galaxies 
might already have an "incorporated mass fraction" compa- 
rable to that of local spiral galaxies (with this fraction likely 
to grow further during subsequent evolution if BX galaxies 
evolve into present-day elliptical galaxies as suggested by 
their clustering; see Adelberge r et al.l2005al) . 

5.3. Dynamical Evolution 

The determination of the ratio of circular velocity v c to z- 
velocity dispersion a may provide an interesting new probe 
of the dynamical evolution of actively star-forming galax- 
ies at high redshift. We consider again the three best cases, 
SSA22a-MD41, Q2343-BX610, and Q2343-BX389, for 
which the modeling of §|4]implies v c /cr and Tin/zi/e of 2-4. 
These ratios are much lower than the value of > 10 for the 
Milky Way disk. We here briefly explore two possible mech- 
anisms through which the local velocity dispersion of the gas 
can be increased in high redshift star-forming galaxies. 

In one scenario, the energy from the gravitational potential 
is converted into random motions during gas infall onto a disk. 
If the conversion is efficient, v c /a gives a measure of the rate 
of gas infall (N. Murray, private communication). Assuming 
that random motions are generated by viscous processes in a 
rotating disk (as considered, e.g ., bv lJog & Ostrikerll9 88 and 
iGammie. Ostriker. & Jogl 19911 for lower velocity dispersions 
in the Galactic disk; see also iGoldreich & Tremaindll978T: 
iFukunagal fl 984 ) and that the rotational and random motion 
energies grow at approximately equal rates, then 

a / Ve nr,<l^, (5) 

V ^acc 



since dispersion is generated on a dynamical timescale and 
rotational energy increases at the mass accretion rate. In 
equation 0, ?dyn = f/v c is the dynamical timescale, f acc = 
M gas /(dM gas /dt) is the mass accretion timescale onto the disk, 
and 77 is a dimensionless parameter with value near unity. In 
this scenario, the values of v c /er - 2-4 for SSA22a-MD41, 
Q2343-BX610, and Q2343-BX389 translate into f acc /f dy n ~ 
4- 16. Given their f dyn ~ 50 Myr at r max ~ 10 kpc, this im- 
plies f acc ~ 200-800 Myr (consistent with scenarios of ac- 
cretion through c old flo ws of baryons from the halo; e.g., 
iDekel & Birnbo im 2006). The rotating disk picture assumed 
here as well as for the modeling of § |4] neglects possible 
contributions to the observed velocity dispersion due to non- 
circular motions, merger-induced flows, or other perturba- 
tions. Therefore, the er/v c values are very likely upper limits 
in the context of equation Q so that the accretion timescale 
estimate s represent almost certainly lower limits. lErb et all 
(2006b) applied population synthesis models to fit the opti- 
cal to NIR (and Spitzer/IRAC mid-IR for a subset) spectral 
energy distribution (SED) of a large BX sample including 
Q2343-BX610 and Q2343-BX389, for which they derived 
stellar ages > 1 Gyr (at 95% confidence), at the high end of 
the r ange of ~ 100 Myr - 1 Gyr typical for BX galaxies (see 
also IShaolev et al.l 12005). Interestingly enough, the lower 
limit for f acc and the stellar ages are consistent with each other. 

An alternative scenario is that of an actively star-forming 
disk in which stellar winds, supernova explosions, and ra- 
diation pressure are responsible for setting the velocity 
dispersion and provid ing vertical support against gravity 
llThompson et al. 12005ft . If the disk reaches margin al stability 
in terms of Toomre's Q-parameter (Toomre 1964), the rela- 
tionship between cr/v c and the gas fraction / gas can be written 
as 

cAc ~ (/gasfi) /a, (6) 

where a = 2 ± 1 is a dimensionless parameter depending on 
the assumed distribution of gas and gravitational potential. In 
this case, the total star formation rate R+ may be described by 
a Schmidt law, 

R± = dM gas /dt = e/ gas M dyn (f dyn ) _1 , (7) 

and the gas exhaustion timescale by star formation is 

h = M gas / (dM m Jdt) = f dyn /e. (8) 

Here e is the star formation efficiency for converting gas into 
stars, which is es timated to be ~ 0. 1 in the Galactic disk (e.g., 
Elme greer]ll999l) . For the inferred range of v c /a ~ 2 — 4, 
equation (jfji and Qrj 1 imply / gas w 0.5— 1. To be as thick as 
inferred from our data the gas disks must be very gas-rich. For 
e 0.1, the gas exhaustion timescale is about 500 Myr. This 
picture, in which the velocity dispersion in a gas-rich star- 
forming disk is related to feedback processes from star for- 
mation, thus leads to estimates of gas exhaustion timescales 
that are plausibly consistent with the stellar ages as well. 

5.4. Disk Stability and Subsequent Evolution 

In view of the above, we now consider the question of the 
lifetime of the present dynamical state of the BX galaxies 
we have observed. In a merger interpretation with two in- 
teracting galaxies within « 10 kpc, the system would form 
a compact (i.e., bulge) merger remnant within w 10 dynami- 
cal timescales, or w 500 Myr. One might suppose that a disk 
configuration would be more stable. However, the disks we 
appear to be observing in some of the BX galaxies differ from 
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local Sa/Sb disk galaxies in that they seem very gas-rich and 
possibly lack a central bulge. Such systems are probably dy- 
n amically unsta b le. 

Ilmmeli et al. I (2004b) present the results of a chemo- 
dynamical simulation of the collapse and fragmentation of 
a gas-rich galactic disk at high redshift (one of the possi- 
ble sc£n^i£sJor diskevolution and bulge formation explored 
by Ilmmeli et al. 1 120043) . In this simulation, Ilmmeli et al- 1 
(2004bj) start from a gaseous disk with a mass inflow rate 
of 120 M yr _1 , growing to a total mass of 1.2 x 10 11 M Q 
in 1 Gyr. This disk fragments into massive star-forming 
clumps, with masses between 0.6 and 7 x 10 9 M Q and (Tga S ~ 
30-50 kms" 1 ~ 0.25 v c . The total mass, the star forma- 
tion rate (reaching a maximum of ~ 220 M Q yr" 1 after ~ 
1 Gyr), the assembly time, and the knotty appearance of this 
disk are very simil ar to the properties of the clu mp-cluster 
galaxies studied by[ Elmegreen & Elmegreeri] ll2005t see also 
Elmeg reen et al.l2004albl) and some of the BX galaxies in our 
sample. The clumpy fragmenting disk is unstable. The star- 
forming clumps sink to the gravitational center by dynamical 
friction against the diffuse underlying disk of gas and form a 
central bulge after ~ 1 Gyr. It is interesting to note that this 
final configuration is similar to those of high redshift galax - 
ies in the hydrodynamical simulations by Abadi et al. (2003), 
and that it also resembles that of a spheroid rather than of a 
disk galaxy. 

In the context of these simulations, it is tempting to spec- 
ulate that some of the BX galaxies in our sample are in the 
mid-stages of the evolution of such unstable gas-rich disks 
later evolving into galaxy bulges. Q2343-BX610 may be in 
a more evolved stage in this sequence, in which an older and 
more metal-rich central bulge of substantial mass has already 
formed. 

5.5. Comparison Between Stellar and Dynamical Masses 

Having estimated dynamical masses for our BX sample, 
an obvious test is to compare t hem w ith photometri c stellar 
mass estimates. Shap levetal.l (120051) and lErb et al.l (|2006b) 
have carried out population synthesis modeling of the SEDs 
of large samples of spectroscopically-confirmed BX galax- 
ies. They found that the averag e stellar mass i n a typical BX 
galaxy is M* « 10 10 3±0 5 M . IShaplev et all (12004) further 
compared the stellar masses and a virial estimate of the dy- 
namical masses of nine K < 20 BX objects, which are at the 
massive tail (M+ ~ 10 11 M Q ) of the population. They con- 
cluded that the stellar masses exceed the dynamical masses by 
a factor of about three for a Salpeter IMF (dN /dm oc m~ ) 
between 0.1 and 100 M Q , which is commonly adopted in 
studies of high redshift galaxies. 

Fro m studies in the Milky Way, it is now clear (e.g.. lScalol 
1986) that the IMF flattens below a few solar masses even 
if a Salpeter slope may be applicable at the high mass end. 
With such a flattening at the l ow mass end, modern versions 
of appropriate local IMFs (e.g. Rroupa 2001; Chabrier 2003) 
lead to total masses ~ 1.6 times lower than the unrealistic 
0. 1 - 100 M Q Salpeter IMF. For a more extreme (but probably 
also unrealistic) 1 - 100 M R Salpeter IMF, total mas ses are 
2. 5 times lo wer. IShaplev et al. I (|2Q04) found that a Kroupa 
or Chabrier IMF is consistent with their virial mass estimates 
derived from the velocity widths for seven BX galaxies. 

The 11 BM/BX SINFONI targets with available NIR pho- 
tometry are included in the sample for which Erb et al. 
(2006b) carried out SED modeling. For a Chabrier IMF, 
their modeling implies and average and rms stellar mass of 



(M*) w io 10 - 8±01 Mq. The dynamical masses we i nferr ed 
from the velocity gradient and line width methods in § !3.2l for 
these objects agree well, with a mean and rms of (Md yn ) rj 
10 io.7±o.2 Mq within t he Ha half-light radius. Taking the re- 
sults of the more detailed kinematic modeling of SSA22a- 
MD41, Q2343-BX610, and Q2343-BX389 as a guide (§|4), 
the ratio of total mass and mass within r\n is about 1.5. 
With this correction, (M dyn )/(M*) is about 10 01±0 ' 3 . For 
Q2343-BX610 and Q2343-BX389, the kinematic model- 
ing gives total M dyn of 9.2 x 10 10 and 2.0 x 10 11 M m , which 
are fact ors of 0.4 and 3 times the M + estimates of E rb et all 
(2006b). These considerations illustrate the uncertainties in- 
volved in deriving photometric stellar masses and dynami- 
cal masses. Nonetheless, the dynamical masses tend to be 
overall ab out equal to the photometric stellar masses with a 
iRroupal or IChabrieJ IMF. Leaving some room for contribu- 
tions from gas and dark matter might imply a slight shortfall 
of the dynamical masses and would then require a still smaller 
fraction of low mass stars or, alternatively, a more top-heavy 
IMF (e.g., as proposed by Baugh et al. 2005). Our dynamical 
mass measurements may thus strengthen the arguments that a 
0.1 - 100 M Q Salpeter IMF is inappropriate for high redshift 
galaxies. 

6. SUMMARY 

We have carried out NIR integral field spectroscopy with 
SINFONI of a sample of 14 UV-selected BM/BX galaxies at 
Z ~ 2 to investigate the morphology and kinematics as traced 
by their Ha line emission. From our Ha linemaps, as well as 
higher resolu tion rest-frame U V imaging from HST for other 
BX galaxies llErb et al. 120 04). it appears that, at least in part, 
the BM/BX galaxies we have observed are not classical ex- 
ponential disks with a central bulge but present rather asym- 
metric and clumpy morphologies. Deep high resolution NIR 
imaging is obviously highly desirable to determine the mor- 
phologies in rest-frame optical continuum emission. 

We detected spatially-resolved velocity gradients in almost 
all galaxies whose integrated size is sufficiently large to allow 
us to make such a measurement. For the six cases for which 
the data provide high quality two-dimensional kinematic in- 
formation, the velocity fields and gradients are smooth and, 
in all but one or two cases, match the expectations for a 
simple rotating disk. The steepest velocity gradients for the 
galaxies in our sample are generally approximately along the 
morphological major axis. For the three best cases among 
our SINFONI sample, we have obtained the first convincing 
asymptotic rotation curves at high redshift, on radial scales 
of > 10 kpc from the dynamical center. The dynamics of the 
BM/BX sample studied here thus appear to be dominated by 
gravitationally-driven motions. 

From these measurements and from kinematic modeling, 
we inferred dynamical masses ranging from ~ 0.5 to 25 x 
10 10 Mq within the Ha half-light radius. Hypothesizing that 
the average v c at a radius of 10 kpc is an appropriate proxy for 
the circular velocity of the dark matter halo implies a virial 
mass for the typical halo of a BM/BX galaxy in our sample 
of 10 1L7±0 - 5 M Q . This is in good agreement with estimates 
based on the cluster ing and co-moving spac e density of the 
BM/BX population (Adelberg er et al. 120053) . The kinematic 
technique may provide a new, independent method for esti- 
mating dark matter halo masses. From the more detailed kine- 
matic modeling of the three best cases, we inferred a ratio of 
peak circular velocity v c to z-velocity dispersion a of ~ 2-4. 
This suggests a mass accretion timescale onto the disk or, al- 
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tentatively, a gas exhaustion timescale of ~ 200 - 800 Myr, 
broadly consistent with the stellar ages. We found that the 
specific angular momenta of the BM/BX objects in our sam- 
ple are similar to those of local late-type galaxies, and that the 
specific angular momenta of their baryonic disks appear to be 
roughly comparable to those of the dark matter halos from 
which they formed. 

We presented arguments, based on the properties of the 
BM/BX galaxies analyzed in this paper, which may suggest 
an interesting scenario for their evolution. Merger-induced 
or (quasi-)adiabatic gas collapse, whereby the specific angu- 
lar momentum of the gas relative to the dark matter halo is 
preserved, has its free energy converted into turbulence in 
the disk, which is one way to explain the low v c /a we ob- 
served. Gas-rich disks growing through infall are likely un- 
stable against fragmentation. The self-gravitating fragments 
result in clumpy and irregular morphologies such as seen in 
the Ha linemaps of several of our BM/BX sources. These 
fragments may transfer angular momentum to infalling gas 
as they sink towards the center of the gravitational potential 
through dynamical friction, leading to the formation of a cen- 
tral bulge. Angular momentum transfer in this way may allow 
the residual gas to keep its high specific angular momentum, 
such that it remains comparable to that of the dark matter halo. 
The collapse of the gas would be relatively inefficient, pre- 



Abadi, M. G., Navarro, J. R, Steinmetz, M., & Eke, V. R. 2003, ApJ, 591, 
499 

Adelberger, K. L., Erb, D. E., Steidel, C. C, Reddy, N. A., Pettini, M., & 

Shapley, A. E. 2005a, ApJ, 620, L75 
Adelberger, K. L., Steidel, C. C, Pettini, M., Shapley, A. E., Reddy, N. A., & 

Erb, D. E. 2005b, ApJ, 619, 697 
Adelberger, K. L., Steidel, C. C, Shapley, A. E., Hunt, M. P., Erb, D. K., 

Reddy, N. A., & Pettini, M. 2004, ApJ, 607, 226 
Armus, L., Heckman, T. M., & Miley, G. K. 1990, ApJ, 364, 471 
Asplund, M., Grevesse, N., Sauval, A. J., Allende Prieto, C, & Kiselman, D. 

2004, A&A, 417, 751 
Baugh, C. M., Lacey, C. G., Frenk, C. S., Granato, G. L., Silva, L., Bressan, 

A., Benson, A. J., & Cole, S. 2005, MNRAS, 356, 1191 
Bonnet, H., et al. 2003, Proc. SPIE, 4839, 329 
Bonnet, H., et al. 2004, The Messenger, 1 17, 17 
Cardone, V. E, & Sereno, M. 2005, A&A, 438, 545 
Chabrier, G. 2003, PASP, 115, 763 

Chapman, S. C, Blain, A. W., Smail, I., & Ivison, R. J. 2005, ApJ, 622, 772 

Cowie, L., Hu, E., & Songaila, A. 1995, AJ, 110, 1576 

Dekel, A. & Birnboim, Y. 2006, MNRAS, in press i astro-ph/04 12300 i 
Dickinson, M., Papovich, C, Ferguson, H. C, & Budavari, T. 2003, ApJ, 
587, 25 

Downes, D., & Solomon, P. M. 1998, ApJ, 507, 615 

Eisenhauer, E, et al. 2003a, Proc. SPIE, 4841, 1548 

Eisenhauer, E, et al. 2003b, The Messenger, 113, 17 

Elmegreen, B.G., 1999, ApJ, 517, 103 

Elmegreen, B. G., Elmegreen, D. M. 2005, ApJ, 627, 632 

Elmegreen, B. G., Elmegreen, D. M., & Hirst, A. C. 2004a, ApJ, 612, 191 

Elmegreen, D. M., Elmegreen, B. G., & Sheets, C. M. 2004b, ApJ, 603, 74 

Elmouttie, M., Koribalski, B., Gordon, S., Taylor, K., Houghton, S., Lavezzi, 

T, Haynes, R., & Jones, K. 1998, MNRAS, 297, 49 
Erb, D. K., Shapley, A. E., Steidel, C. C, Pettini, M., Adelberger, K. L., Hunt, 

M. P., Moorwood, A. F. M., & Cuby, J.-G. 2003, ApJ, 591, 101 
Erb, D. K., Steidel, C. C, Shapley, A. E., Pettini, M., & Adelberger, K. L. 

2004, ApJ, 612, 122 
Erb, D. K., Shapley, A. E., Pettini, M., Steidel, C. C, Re ddy, N. A., & 

Adelberger, K. L. 2006a, ApJ, in press i astro-ph/0602473 i 
Erb, D. K., Steidel, C. C, Shapley, A. E., Pettini, M., Reddy, N. A., & 

Adelberger, K. L. 2006b, ApJ, submitted 
Fan, X., et al. 2001, AJ, 121, 54 
Forster Schreiber, N. M., et al. , 2004, ApJ, 616, 40 
Fontana, A., et al. 2003, ApJ, 594, L9 
Franx, M., et al. , 2003, ApJ, 587, L79 
Fukunaga, M. 1984, PASP, 36, 433 

Gammie, C. F., Ostriker, J. P., & Jog, C. J. 1991, ApJ, 378, 565 



serving the relatively low fraction of baryonic mass that we 
inferred for our sample. Although some of these arguments 
are speculative, they may suggest a plausible and consistent 
picture that may explain a number of features observed for 
some of the BM/BX galaxies studied here. 

The results presented in this work apply to the BM/BX sam- 
ple that we have observed with SINFONI. Admittedly, the 
sample may not be representative of the BM/BX population as 
a whole because the targets were preferentially selected to be 
large, or to exhibit significant velocity gradients or large line 
widths based on existing long-slit spectroscopy. Nonetheless, 
the SINFONI data we have obtained show the improvement 
in our ability to investigate the morphologies and kinematics 
of high redshift galaxies afforded by NIR integral field spec- 
troscopy. It will undoubtedly be very interesting to extend 
such studies to other high redshift populations and further in- 
vestigate the issues addressed in this work in the future. 
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TABLE 1 

Summary of the SINFONI observations of the BM/BX sample 



Source 


ZH Q a 


Band 


Scale 


t- b 
'nit 

(s) 


Run c 


Q 1 307 — rJM 1 1 oi 


1.4101 


H 


0. 125 


14400 


Maru5 


Qloij — BXi/oa 


2.4087 


K 


Of 125 


15600 


Mar05, Apr05 


Q1623 — BX455 


2.4071 


K 


Of 125 


12000 


Mar05 


Q1623 — BX502 


2.1555 


K 


Of 05 


9600 


Apr05 


y lOZJ — DAJZO 


2 2684 






1 JJW 


Jul04 


Q1623-BX599 


2.3313 


K 


Of 125 


5400 


Jul04 


Q1623-BX663 


2.4335 


K 


0"125 


9000 


Jul04 


SSA22a-MD41 


2.1710 


K 


0"125 


28800 


Nov04, Jun05 


Q2343-BX389 


2.1728 


K 


0"125 


18000 


Oct05 


Q2343-BX610 


2.2102 


K 


0"125 


10800 


Jun05 






H 


Of 125 


18000 


Aug05, Oct05 






.1 


0"125 


14400 


Oct05 


Q2346-BX404 


2.0298 


K 


0"125 


8100 


Jul04 


Q2346-BX405 


2.0308 


K 


0"125 


8100 


Jul04 


Q2346-BX416 


2.2406 


K 


0"125 


7200 


Dec04 


Q2346-BX482 


2.2569 


K 


Of 125 


7200 


Nov04 



a Systemic vacuum redshift derived from the integrated Ha line emission in the SINFONI data (as 
described in §[5}. Q2348-BX404 and BX405 form a pair, with angular separation of 3." 6 (or 30 kpc 
at their redshift). 

b Total on-source integration time. 

c SINFONI observing runs where the data were taken: Jul04: commissioning run on 2004 luly 9 to 
18; Nov04: GTO run on 2004 November 29 and 30; Dec04: GTO run on 2004 December 20 to 22; 
Mar05: GTO run on 2005 March 13 to 22; Apr05: GTO run on 2005 April 4 to 7. Jun05: GTO run 
on 2005 lune 15 to 17. Aug05: GTO run on 2005 August 27 to September 4. Oct05 : GTO run on 
2005 October 2 to 12. 



TABLE 2 

Properties derived from the SINFONI observations of the BM/BX sample. 



Source 


ZHo' 


[N II]/Ha b 


n/2 c 


Av d 




sin 2 (i) X M dyn e 








(kpc) 


(kms-') 


(kms-') 


(10 10 M Q ) 


Q1307-BM1163 


1.4101 


0.22 ± 0.03 


3.1 


310 


60 


0.84 


Q1623-BX376a 


2.4087 


0.15±0.1 


3.1 


210 


30 


0.33 


Q1623-BX455 


2.4071 


0.14 ±0.07 


3.5 


270 


55 


0.73 


Q1623-BX502 


2.1555 


<0.05 


2.0 


170 


50 


0.21 


Q1623-BX528 


2.2684 


0.32 ±0.05 


6.6 


280 


65 


1.7 


Q1623-BX599 


2.3313 


0.24 ±0.08 


3.4 


375 


50 


1.1 


Q1623-BX663 


2.4335 


0.24 ±0.08 


5.3 


360 


95 


2.4 


SSA22a-MD41 


2.1710 


0.15 ±0.05 


6.0 


300 


140 


3.2 


Q2343-BX389 


2.1728 


0.23 ± 0.07 


7.0 


570 


235 


11 


Q2343-BX610 


2.2102 


0.41 ±0.02 


5.4 


385 


165 


4.6 


Q2346-BX404 


2.0298 


0.20 ±0.06 


2.8 


210 


20 


0.27 


Q2346-BX405 


2.0308 


<0.05 


4.6 


170 


30 


0.37 


Q2346-BX416 


2.2406 


0.17 ±0.05 


2.8 


315 


70 


0.84 


Q2346-BX482 


2.2569 


0.19 ±0.05 


5.8 


300 


100 


2.3 



a Systemic vacuum redshift derived from the integrated Ha line emission. 

b Ratio of the integrated [N II] A6584 and Ha line fluxes. For Q1307-BM1163, Q1623-BX599, 
and Q1623 -BX663, there is a broad component to Ha, making estimates of the [N II]/Ha ratio 
more uncertain. 

c Rest-frame intrinsic half-light radius n/2 along the major axis of the detected Ha light distribution 
(see §|5J- The observed HWHM along the morphological major axis of each galaxy was corrected 
for a typical NIR seeing with FWHM = Of 5 except for Q1623-BX502, for which AO-assisted 
observations were taken, resulting in an effective resolution with FWHM = Of 3. 

d Intrinsic FWHM velocity width Av and half-total velocity spread v r = (v max - Vmin)/2 over the 
Ha emitting regions (see § The FWHM line widths for each galaxy were corrected for the 
instrumental resolution of ~ 80 kms~' . 

c Inferred dynamical mass within n/2 (uncorrected for inclination), averaged from the "line width" 
and "velocity gradient" methods, as described in §|5| 
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Fig. 1. — Ha morphology and kinematics of the sample of 14 BM/BX objects observed with SINFONI. For each galaxy, we show the 
velocity-integrated Ha linemap (top insets) and a position-velocity (p-v) diagram extracted in a w 0."5-wide synthetic slit (bottom insets). 
The linemaps are displayed with a linear intensity scale increasing from black to white, and the p-v diagrams are colour-coded with a linear 
intensity scale increasing from magenta to red. Each linemap is labeled with the source name and redshift, and the position of the synthetic slit 
is shown as well (thin dotted line). In the inset pair at the bottom right, the angular and velocity scales for the linemaps and position-velocity 
diagrams are marked, and the full synthetic slit width used for all objects is shown. The angular size of 1" is also indicated in all linemaps 
(thick solid line), corresponding to ~ 8.3 kpc at the redshifts of the objects. For all sources, the data were obtained with the 0"l25 pixel -1 
scale, with the exception of Q1623-BX502, which was observed with the 0"05 pixel -1 scale with AO correction on a nearby star. In the 
p-v diagrams, one pixel corresponds to 0.000195 /im along the horizontal (wavelength) axis for the /f-band data of Q1307-BM1163 and to 
0.000245 fim for all other sources, observed in the A"-band. In all linemaps, north is up and east is to the left. The objects show a variety 
of Ha morphologies and velocity gradients. The [N II] A 6584 emission line is visible and labeled in the p-v diagrams of Q1307-BM1163, 
Q2343-BX389, Q2343-BX610, Q2346-BX416, and Q2346-BX404. 



Integral field spectroscopy of z ~ 2 galaxies 



15 



QI623-BX663 




SSA22-MD41 






Fig. 2. — Two-dimensional Ha kinematics for five of our BM/BX sample galaxies observed with SINFONI. Each row corresponds to the 
results for one source. The panels to the left show the integrated Ha linemap, with the angular scale of 1" indicated (vertical bar), corresponding 
to « 8.3 kpc at the redshifts of the sources. The middle panels show the spatial distribution of velocity v (relative to the systemic velocity), 
obtained from Gaussian fits to the data cubes after smoothing spatially with a two-dimensional Gaussian of FWHM = 3 pixels (or 0"38, resulting 
in an effective resolution of ~ 0"6). The panels to the right show the spatial distribution of the FWHM velocity width obtained from the same 
fits. The velocity and FWHM maps are color-coded with a linear stretch and such that values increase from magenta to red; minimum and 
maximum values for each galaxy are labeled on the maps. An outline of the integrated Ha linemap is overplotted on the velocity and FWHM 
maps. In all images, north is up and east is to the left. 
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Fig. 3. — Rest-frame optical emission and Ha kinematics of Q2343-BX610. The maps shown are the integrated Ha line emission {top left), 
the integrated [N II] A6584 line emission {top center), the integrated line-free observed K continuum emission {top right), the [N II] A6584/Ha 
line ratio {bottom left), the Ha velocity field v(Ha) {bottom center), and the Ha FWHM velocity width Av(Ha) {bottom right). The vertical 
bar left of the Ha linemap is 1" in length, or 8.26 kpc at the redshift of Q2343-BX610. Each map is color-coded with a linear scaling and 
such that the values increase from magenta to red. The minimum and maximum values of [N II] A6584/Ha ratio (unitless), and of v(Ha) and 
Av(Ha) (in kms~') are labeled below their respective maps. All maps have been extracted from the reduced data cube after smoothing spatially 
by a two-dimensional Gaussian of FWHM = 3 pixels (or 0."38). Contours of the Ha line emission are overplotted on each map. In all images, 
north is up and east is to the left. 
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Fig. 4.— Rotation curves of SSA22a-MD41, Q2343-BX610, and Q2343-BX389 derived from the Ha line emission. In each panel, 
the data points show the rotation velocities (relative to the systemic velocity) as a function of position along the kinematic major axis of each 
galaxy. The velocities were obtained from Gaussian fits to the line profiles measured in spectra extracted in synthetic circular apertures. The 
horizontal error bars indicate the diameter of the synthetic apertures. The vertical error bars correspond to the formal la uncertainties of the 
velocities. The rotation curves of the disk models that best fit the observed Ha kinematics of each galaxy and described in §[4]are also plotted. 
The best-fit ring models are plotted as red solid lines, and the best-fit exponential disk models (with central hole) are shown as black dashed 
lines. 




Fig. 5. — Two-dimensional Ha velocity field of Q2343-BX610. The velocity field derived from the observed Ha line emission is shown in 
colours, with a linear scaling increasing from magenta to red. The superpo sed c ontours show the isovelocity map from the best-fit rotating disk 
model with values relative to the systemic velocity in kms~' labeled (see § 14.31 . 



